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!>■ . Abstract 

^ ' The structure of events in high-energy collisions is complex and not predictable 



oo 



from first principles. Event generators allow the problem to be subdivided into 
more manageable pieces, some of which can be described from first principles, 
while others need to be based on appropriate models with parameters tuned to 
data. In these lectures we provide an overview, discuss how matrix elements 
■ are used, introduce the machinery for initial- and final-state parton showers, 

! explain how matrix elements and parton showers can be combined for optimal 

CN I accuracy, introduce the concept of multiple parton-parton interactions, com- 

ment briefly on the hadronization issue, and provide an outlook for the future. 

^ ' 

^ . 1 Introduction 

p • Given the current landscape in experimental high-energy physics, these lectures are focused on appli- 

cations of event generators for hadron colliders like the Tevatron and LHC. Much of it would also be 
D I relevant for e+e" machines like LEP and ILC or e^p machines like HERA, but with some differences 

' not discussed here. Heavy-ion physics is not at all addressed, since it involves rather different aspects, 

specifically the potential formation of a quark-gluon plasma. Further, within the field of high-energy 
^ ' PP/PP collisions, the emphasis will be on the common aspects of QCD physics that occurs in all col- 

^ I lisions, rather on those aspect that are specific to a particular physics topic, such as B production or 

supersymmetry. Both heavy ions and other physics topics are instead covered by other lectures at this 
school. 

Section 2 contains a first overview of the physics picture and the generator landscape. Thereafter 
section 3 describes the usage of matrix elements, section 4 the important topics of initial- and final-state 
showers, and section 5 how showers can be matched to different hard processes. The issue of multiple 
interactions and their role in mimimum-bias and underlying-event physics is introduced in section 6, 
followed by some brief comments on hadronization in section 7. The article concludes with an outlook 
on the ongoing generator-development work in section 8. 

Slides for these and other similar lectures [1] are complementary to this writeup in style and 
contents, including many (colour) illustrations absent here. Other useful resources include the "Les 
Houches Guidebook to Monte Carlo Generators for Hadron Collider Physics" [2] and a recent review on 
QCD physics at the Tevatron and LHC [3]. 

2 Overview 

In real life, machines produce events that are stored by the data acquisition system of a detector. In 
the virtual reality, event generators like Herwig [4] and Pythia [5] play the role of machines like the 
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Tevatron and LHC, and detector simulation programs like Geant 4 [6] the role of detectors like ATLAS 
or CMS. The real and virtual worlds can share the same event reconstruction framework and subsequent 
physics analysis. It is by understanding how an original physics input is distorted step-by-step in the 
better-controlled virtual world that an understanding can be gained of what may be going on in the real 
world. For approximate studies the detector simulation and reconstruction steps can be shortcut, so that 
generators can be used directly in the physics studies. 

A number of physics analyses would not be feasible without generators. Specifically, a proper 
understanding of the (potential) signal and background processes is important to separate the two. The 
key aspect of generators here is that they provide a detailed description of the final state so that, ideally, 
any experimental observable or combination of observables can be predicted and compared with data. 
Thereby generators can be used at various stages of an experiment: when optimizing the detector and 
its trigger design to the intended physics program, when estimating the feasibility of a specific physics 
study, when devising analysis strategies, when evaluating acceptance corrections, and so on. 

However, it should always be kept in mind that generators are not perfect. They suffer from 
having to describe a broad range of physics, some of which is known from first principles, while other 
parts are modelled in different frameworks. (In the latter case, a generator actually acts as a vehicle of 
ideology, where ideas are disseminated in prepackaged form from theorists to experimentalists.) Given 
the limited resources, different authors may also have invested more or less time on specific physics 
topics, and therefore these may be more or less well modelled. It always pays to shop around, and to 
compare several approaches before drawing too definite conclusions. BUnd usage of a generator is not 
encouraged: then you are the slave rather than the master. 

Why then Monte Carlo event generators? Basically because Einstein was wrong: God does throw 
dice! In quantum meachanics, calculations provide the probability for different outcomes of a measure- 
ment. Event-by-event, it is impossible to know beforehand what will happen: anything that is at all 
allowed could be next. It is only when averaging over large event samples that the expected probability 
distributions emerge — provided we did the right calculation to high enough accuracy. In generators, 
(pseuo)random numbers are used to make choices intended to reproduce the quantum mechanical prob- 
abilities for different outcomes at various stages of the process. 

The buildup of the structure in an event occurs in several steps, and can be summarized as follows: 

- Initially two hadrons are coming in on a collision course. Each hadron can be viewed as a bag of 
partons — quarks and gluons. 

- A collision between two partons, one from each side, gives the hard process of interest, be it for 
physics within or beyond the standard model: ug ug, ud , gg h^, etc. (Actually, 
the bulk of the cross section results in rather mundane events, with at most rather soft jets, or of 
a simple elastic or diffractive character that is not easily described as partonic processes. Such 
events usually are filtered away at an early stage, however.) 

- When short-lived "resonances" are produced in the hard process, such as the top, or Z^, their 
decay has to be viewed as part of this process itself, since e.g. spin correlations are transferred 
from the production to the decay stages. 

- A collision impUes accelerated colour (and often electromagnetic) charges, and thereby brems- 
strahlung can occur. Emissions that can be associated with the two incoming colHding partons are 
called Initial-State Radiation (ISR). As we shall see, such emissions can be modelled by so-called 
space-like parton showers. 

- Emissions that can be associated with outgoing partons are instead called Final-State Radiation 
(FSR), and can be approximated be time-Uke parton showers. Often the distinction between a hard 
process and ISR and FSR is ambiguous, as we shall see. 

- So far we only extracted one parton from each incoming hadron to undergo a hard collision. But 
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the hadron is made up of a multitude of further partons, and so further parton pairs may collide 
within one single hadron-hadron collision — multiple interactions (MI). (Not to be confused with 
pileup events, when several hadron pairs colUde during a bunch-bunch crossing, but with obvious 
analogies.) 

- Each of these further collisions also may be associated with its ISR and FSR. 

- The colliding partons take a fraction of the energy of the incoming hadrons, but much of the energy 
remains in the beam remnants, which continue to travel essentially in the original directions. These 
remnants also carry colours that compensate the colour taken away by the colliding partons. 

- At short timescales, when partons are close to each other, the principle of asymptotic freedom 
tells us that we can think of each parton as freely moving along its trajectory. However, as time 
goes by and the partons recede from each other, confinement forces become significant. The struc- 
ture and time evolution of these force fields cannot be described from first principles within any 
calculational technique currently at our disposal, so models have to be introduced. One common 
approach is to assume that a separate confinement field is stretched between each colour and its 
matching anticolour, with each gluon considered as a simple sum of a colour and an anticolour, 
and all colours distinguishable from each other (the Nc oo limit). 

- Such fields can break up by the production of new quark-antiquark pairs that screen the endpoint 
colours, and where a quark from one break (or from an endpoint) can combine with an antiquark 
from an adjacent break to produce a primary hadron. This process is called hadronization. 

- Many of those primary hadrons are unstable and decay further at various timescales. Some are 
sufficiently long-hved that their decays are visible in a detector, or are (almost) stable. Thereby 
we have reached scales where the event-generator description has to be matched to a detector- 
simulation framework. 

- It is only at this stage that experimental information can be obtained and used to reconstruct back 
what may have happened at the core of the process, e.g. whether a Higgs particle was produced or 
not. 

The Monte Carlo method allows these steps to be considered sequentially, and within each step 
to define a set of rules that can be used iteratively to construct a more and more complex state, maybe 
ending with hundreds of particles moving out in different directions. Since each particle contains of the 
order of ten degrees of freedom (flavour, mass, momentum, production vertex, lifetime, . . . ) we realize 
that thousands of choices are involved for a typical event. The aim is to have a sufficiently reahstic 
description of these choices that both the average behaviour and the fluctuations around this average are 
well decribed. 

Schematically, the cross section for a range of final states is provided by 

'-'^final state '-'^hard process ^tot,hard process— >final state ; (1) 

properly integrated over the relevant phase-space regions and summed over possible "paths" (of show- 
ering, hadronization, etc.) that lead from a hard process to the final state. That is, the dimensional 
quantities are associated with the hard process; subsequent steps are handled in a probabiUstic approach. 

The spectrum of event generators is very broad, from general-purpose ones to more specialized 
ones. Herwig and Pythia are the two most commonly used among the former ones, with ISAJET [7] 
and Sherpa [8] as the other two main programs in this category. Among more specialized programs, 
many deal with the matrix elements for some specific set of processes, a few with topics such as parton 
showers or particle decays, but there are e.g. no freestanding programs that handle hadronization. In the 
end, many of the specialized programs are therefore used as "plugins" for the general-purpose ones. 



3 Matrix elements and their usage 

From the Lagrangian of a theory the Feynman rules can be derived, and from them matrix elements 
can be calculated. Combined with phase space it allows the calculation of cross sections. As a simple 
example consider the scattering of quarks in QCD, say u{l) d{2) n(3) d(4), a process similar to 
Rutherford scattering but with gluon exchange instead of photon ditto. The Mandelstam variables are 
defined as s = {pi + V'l?'^ i = {Pi — Ps)^ and u = {pi — Pa)'^. In the cm frame of the colUsion s is the 
squared total energy and i,u = —s{l ^ cos 9)/2 where 9 is the scattering angle. The differential cross 
section is then 

da TT 4 2 + 

which diverges roughly Uke dp'^/p\ for transverse momentum p± —>■ 0. We will come back to this issue 
when discussing multiple interactions; for now suffice to say that some lower cutoff p_Lmin need to be 
introduced. Similar cross sections, differing mainly by colour factors, are obtained for g 51 qg and 
9 9 ~^ 9 9- few further QCD graphs, like g g qq, are less singular and give smaller contributions. 
These cross sections then have to be convoluted with the flux of the incoming partons i and j in the two 
incoming hadrons A and B: 
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diidx2d«7i'"(xi,Q=)/f' (X2,(J2) ^ . (3) 
at 

The parton density functions (PDFs) of gluons and sea quarks are strongly peaked at small momentum 
fractions xi Ei/EA, X2 ~ Ej/Eb- This further enhances the peaking of the cross section at small p± 
values. Nevertheless, with high machine luminosity the jet cross section can be studied out to quite high 
values. 

The cross section of other processes can be suppressed by two main effects. Firstly, for massive 
particles the p± spectrum is strongly dampened below the respective mass scale, and it is only above it 
that these precesses have a chance to stand up above the QCD background. Secondly, the processes may 
involve electroweak (or other small) couplings rather than strong ones. 

In order to address the physics of interest a large number of processes, both within the Standard 
Model and in various extensions of it, have to be available in generators. Indeed many can also be found 
in the general-purpose ones, but by far not enough. Further, often processes are there available only to 
lowest order, while experimental interest may be in higher orders, with more jets in the final state, either 
as a signal or as a potential background. So a wide spectrum of matrix-element-centered programs are 
available [9], some quite specialized and others more generic. 

The way these programs can be combined with a general-purpose generator is illustrated in Fig. 1. 
In the study of Supersymmetry (SUSY) it is customary to define a model in terms of a handful parameters, 
e.g. specified at some large Grand Unification scale. It is then the task of a spectrum calculator to turn 
this into a set of masses, mixings and couplings for the physical states to be searched for. Separately, 
the matrix elements can be calculated with these properties as unknown parameters, and only when 
the two are combined is it possible to speak of physically relevant matrix-element expressions. These 
matrix elements now need to be combined with PDF's and sampled in phase space, preferable with 
some preweighting procedure so that regions of phase space with high cross sections are sampled more 
frequently. The primarily-produced SUSY particles typically are unstable and undergo sequential decays 
down to a lightest supersymmetric particle (LSP), again with branching ratios and angular distributions 
that should be properly modelled. The LSP would be neutral and escape undetected, while other decay 
products would be normal quarks and leptons. 

It is at this stage that general-purpose programs take over. They describe the showering associ- 
ated with the above process, the presence of additional interactions in the same hadron-hadron collision, 
the structure of beam remnants, and the hadronization and decays. They would still rely on the exter- 
nally suppUed PDF's, and potentially make use of programs dedicated to r and B decays, where spin 
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Fig. 1: Example how different programs can be combined in the event-generation chain. 



information and form factors require special encoding. Even after the event has been handed on to the 
detector-simulation program some parts of the generator may be used in the simulation of secondary 
interactions and decays. 

Several standards have been developed to further this interoperability. The Les Houches Accord 
(LHA) for user processes [10] specifies how parton-level information about the hard process and sequen- 
tial decays can be encoded and passed on to a general-purpose generator. Originally it was defined in 
terms of two Fortran commonblocks, but more recently a standard Les Houches Event File format [1 1] 
offfers a language-independent alternative approach. The Les Houches Accord Parton Density Functions 
(LHAPDF) library [12] makes different PDF sets available in a uniform framework. The SUSY Les 
Houches Accord (SLHA) [13] allows a standardized transfer of masses, mixings, couplings and branch- 
ing ratios from spectrum calculators to other programs. Finally, the HepMC C-i-i- event record [14] 
succeeds the HEPEVT Fortran one [15] as a standard way to transfer information from a generator on to 
the detector-simulation stage. One of the key building blocks for several of these standards is the PDG 
codes for all the most common particles [16], also in some scenarios for physics beyond the Standard 
Model. 

The 2^2 processes we started out with above are about the simplest one can imagine at a hadron 
collider. In reality one needs to go on to higher orders. lnO{a^) two new kind of graphs enter. One kind 
is where one additional parton is present in the final state, i.e. 2 — > 3 processes. The cross section for 
such processes is almost always divergent when one of the parton energies vanish (soft singularities) or 
two partons become collinear (collinear singularities). The other kind is loop graphs, with an additional 
intermediate parton not present in the final state, i.e. a correction to the 2 — > 2 processes. Strictly 
speaking, at 0(ag) one picks up the interference between the lowest-order graph and the loop graph, and 
this interference has negative divergences that exactly cancel the positive ones above, with only finite 
terms surviving. For inclusive event properties such next-to-leading order (NLO) calculations lead to 
an improved accuracy of predictions, but for more exclusive studies the mathematical cancellation of 
singularities has to be supplemented by more physical techniques, which is far from trivial. 

The tricky part of the calculations is the virtual corrections. NLO is now state-of-the-art, with 
NNLO still in its infancy. If one is content with Bom-level diagrams only, i.e. without any loops, it 
is possible to go to quite high orders, with up to something like eight partons in the final state. These 
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Fig. 2: The "factorization" of a 2 ^ n process. 

partons have to be kept well separated to avoid the phase-space regions where the divergences become 
troublesome. In order to cover also regions where partons become soft/coUiner we therefore next turn 
our attention to parton showers. 

4 Parton showers 

To iterate, the emission rate for a branching such as g — > qg diverges when the gluon either becomes 
coUinear with the quark or when the gluon energy vanishes. The QCD pattern is similar to that for e — > 67 
in QED, except with a larger coupling, and actually a coupling that also increases for smaller relative 
p± in a branching, thereby further enhancing the divergence. Furthermore the non-Abelian character of 
QCD leads to g ^ gg branchings with similar divergences, without any correspondence in QED. The 
third main branching, g ^ qq with its 7 e^e~ QED equivalence, does not have the soft divergence 
and is less important. 

Now, if the rate for one emission of a gluon is big, then also the rate for two or more will be big, 
and thus the need for high orders and many loops in matrix-element-based descriptions. With showers 
we introduce two new concepts that make like easier, 

(1) an iterative structure that allows simple expressions for q ^ qg, g ^ gg and g ^ qq branchings to 
be combined to build up complex multiparton final states, and 

(2) a Sudakov factor that offers a physical way to handle the cancellation between real and virtual diver- 
gences. 

Neither of the simplifications is exact, but together they allow us to provide sensible approximate answers 
for the structure of emissions in soft and collinear regions of phase space. 

4.1 The shower approach 

The starting point is to "factorize" a complex 2 n process, where n represents a large number of 
partons in the final state, into a simple core process, maybe 2 — > 2, convoluted with showers. Fig. 2. 
To begin with, in a simple ud — > ud process the incoming and outgoing quarks must be on the mass 
shell, i.e. satisfy = i?^ — = ~ 0, at long timescales. By the uncertainty principle, however, 
the closer one comes to the hard interaction, i.e the shorter the timescales considered, the more off-shell 
the partons may be. (In the oldfashioned perturbative language particles were always on mass-shell, and 
the uncertainty relation allowed energy not to be conserved temporarily. In the modern Feynman-graph 
language four-momentum is conserved at each vertex, but intermediate "propagator" particles need not 
be on the mass shell. The final physics is the same in both languages.) 



6 



Thus the incoming quarks may radiate a succession of harder and harder gluons, while the out- 
going ones radiate softer and softer gluons. One definition of hardness is how off-shell the quarks are, 
_ _ p2|^ but we will encounter other variants later. In the initial-state radiation (ISR) part of 

the cascade these virtualities are spacelike, < 0, hence the alternative name spacelike showers. Cor- 
respondingly the final-state radiation (FSR) is characterized by timelike virtualities, > 0, and hence 
also called timeUke showers. 

To see where this distinction comes from, consider the kinematics of an arbitrary branching a 
be, with a defined to be moving along the +z axis. Then it is useful to introduce Ughtcone momenta 
p± = E ± pz, so that the relation p^ = E'^ — p^ — p^ — pi = m? translates to p+p- = m? + + Py = 
m? +p]_. Now define the splitting of p+ by p+b = zp+a and p+c = (1 — z)p+a- Obviously p_Lc = — p±6 
so p\ = p\fj = p^^. Remains to ensure conservation of p- = {vn? +p\)/p+'- 

P-a=P-b+P-c ^ — = ' + n ^ ^ ^a = — + r^+ n X ■ (4) 

p+a Zp+a {l-z)p+a Z l-Z Z[l - Z) 

In an initial-state branching the incoming a should be (essentially) massless, and if c does not interact 
any further it should also be massless. This gives mf = —{1 — z)p\ < 0, a virtuaUty that is acceptable 
if h is on its way in to a hard scattering, i.e. is destined only to live for a short while. In ISR the Qf 
virtuahties, such as Q\ and in Fig. 2, are usually defined as — m| to keep them positive definite. For 
a final-state branching, assume that h and c will not branch any further and thus are massless, while a is 
an intermediate particle coming from the hard interaction. Then m^ = p\/ {z{l — z)) > , cf. and 
Q\ in Fig. 2, here with Qf = +m^. 

The cross section for the whole 2 ^ n graph is associated with the cross section of the hard 
subprocess, with the approximation that the other virtuahties can be neglected in the matrix-element 
expression. In the hmit that all the <C this should be a good approximation. In other words, first 
the hard process can be picked without any reference to showers, and only thereafter are showers added 
with unit probability. But, of course, the showers do modify the event shape, so at the end of the day the 
cross section is affected. For instance, the total transverse energy £^_Ltot of an event is increased by ISR, 
so the cross sections of events with a given £^_Ltot is increased by the influx of events that started out with 
a lower E'j.tot in the hard process. 

It is important that the hard-process scale is picked to be the largest one, i.e. > Qf in 
Fig. 2. If e.g. Q\ > then instead the ug ug subgraph ought to be chosen as hard process, and the 
gluon of virtuality ought to be part of the ISR off the incoming d. Without such a criterion one might 
doublecount a given graph, or even count it once for every possible subgraph inside the complete 2 ^ n 
graph. In addition, the approximation of neglecting virtualities in the hard-scattering matrix elements 
obviously becomes worse the more the incoming and outgoing partons are off-shell, another reason not 
to put a larger scale than necessary in the shower part. 

4.2 Final-state radiation 

Let us next turn to a more detailed presentation of the showering approach, and begin with the simpler 

final-state stage. This is most cleanly studied in the process e~^e~ — > 7*/Z° qq. The first-order 
correction here corresponds to the emission of one additional gluon, by either of the two Feynman graphs 
in Fig. 3. Neglect quark masses and introduce energy fractions xj = 2Ej/Ecra in the rest frame of the 
process. Then the cross section is of the form 

dcTME "s 4 xl + xl 

(To 27r 3 (1 — — a;2) 
where ctq is the qq cross section, i.e. without the gluon emission. 

Now study the kinematics in the limit X2 — > 1. Since 1 — X2 = m^^/E^^ we see that this 
corresponds to the "colhnear region", where the separation between the q and g vanishes. Equivalently, 
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the virtuality = Q\ = of the intermediate quark propagator i in Fig. 3a vanishes. Although the 
full answer contains contributions from both graphs it is obvious that, in this region, the amplitude of 
the one in Fig. 3a dominates over the one in Fig. 3b. We can therefore view the process as j*/Z^ — > qq 
followed by g ^ qg. Define the energy sharing in the latter branching by Eq = zEi and Eg = {1 — z)Ei. 
The kinematics relations then are 

l-X2 = -^ = -^ ^ dX2 = -^ (6) 

-^cm -'--'cm -^cm 

xi ^ z dxi dz (7) 

X3 ~ 1 - 2 (8) 

so that 

_ ddME _ Os dx2 4 xl + xj _ Os dQ2 4 1 + z^ 
' ~ ao ~ 2n (1 - X2) 3 (1 - xi) "^^^ ^ 27r 3 1-z ^ 

Here dQ^/Q^ corresponds to the "collinear" or "mass" singularity and dz/{l — z) = dEg/Eg to the 
soft-gluon singularity. 

The interesting aspect of eq. (9) is that it is universal: whenever there is a massless quark in the final 
state, this equation provides the probability for the same final state except for the quark being replaced 
by an almost collinear qg pair (plus some other slight kinematics adjustments to conserve overall energy 
and momentum). That is reasonable: in a general process any number of distinct Feynman graphs may 
contribute and interfere in a nontrivial manner, but once we go to a collinear region only one specific 
graph will contribute, and that graphs always has the same structure, in this case with an intermediste 
quark propagator. Corresponding rules can be derived for what happens when a gluon is replaced by a 
collinear gg or qq pair. These rules are summarized by the DGLAP equations [17] 



27r Q2 



dVa^bc = ^^Pa^bc{z)dz (10) 



4 1 + ^2 
3 I- z ' 
(1 -z{l-z )f 
z{l-z) 



where Pq->qg = tt ' (H) 

^g^gg = 3 ^^ ,Z . (12) 



ft f 

Pg-tqq = ^ (^^ + (1 — z)'^) {uf = HO. of quark flavours) . (13) 

Furthermore, the rules can be combined to allow for the successive emission in several steps, e.g. where 
a. q ^ qg branching is followed by further branchings of the daughters. That way a whole shower 
develops. Fig. 4. 

Such a picture should be reliable in cases where the emissions are strongly ordered, i.e. » 
Q2 ^ Qs- ■ ■■ Showers would not be useful if they only could be applied to strongly-ordered parton 
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Fig. 4: A cascade of successive branchings. 



configurations, however. A further study of the ^* /Z^ qqg example shows that the simple sum of 
the q qg and q ^ qg branchings reproduce the full matrix elements, with interference included, to 
better than a factor of 2 over the full phase space. This is one of the simpler cases, and of course one 
should expect the accuracy to be worse for more complicated final states. Nevertheless, it is meaningful 
to use the shower over the whole strictly-ordered, but not necessarily strongly-ordered, region Qi > 
Q2 > . . .to obtain an approximate answer for multiparton topologies for which the complete matrix 
elements would be too lengthy. 

With the parton-shower approach, the big probability for one branching q ^ qg turns into a big 
probability for several successive branchings. Nevertheless we did not tame the fact that probabilities 
blow up in the soft and coUinear regions. For sure, perturbation theory will cease to be meaningful at 
so small scales that as{Q^) diverges; there confimenent effects and hadronization phenomena take 
over. Typically therefore some lower cutoff at around 1 GeV is used to regulate both soft and coUinear 
divergences: below such a scale no further branchings are simulated. Whatever perturbative effects may 
remain are effectively pushed into the parameters of the nonperturbative framework. That way we avoid 
the singularities, but we can still have "probabilities" well above unity, which does not seem to make 
sense. 

This brings us to the second big concept of this section, the Sudakov (form) factor [18]. In the 
context of particle physics it has a specific meaning related to the properties of the loop diagrams, but 
more generally we can just see it as a consequence of the conservation of total probability 



P (nothing happens) = 1 — P (something happens) 



(14) 



where the former is multiplicative in a time-evolution sense: 



7'nothing(0 <t<T)= Pnothing(0 < t < Ti) Nothing (Tl < t < T) . 

Now subdivide further, with Tj = {i/n)T, < i < n: 



(15) 



n-1 



Nothing (0 <t<T) 



lim TT Pnothing(Ti < t < Tj+i) 



j=0 



n-1 



lim TT (1 



^something (^i ^ t ^ ^i+l)) 



i=0 




q 



=^ dV^r.t{T) = dnomething(T) exp(-^^^^^^^l^^dt) . (16) 

That is, the probability for something to happen for the first time at time T is the naive probability 
for this to happen, times the probability that this did not yet happen. As such it applies to a host of 
situations. Take the example of football (relevant at the time of the school). Assume that players are 
equally energetic and skillful from the first minute of the match to the last. Then the chances of scoring a 
goal is uniform in time, but the probability of scoring the^r^? goal of the match is bigger at the beginning, 
because later on any goal could well be the second or third. 

In physics a common example is that of radioactive decay. If the number of undecayed radioactive 
nuclei at time t is J\f{t), with initial number A/q at time t = 0, then a naive ansatz would be dM /dt = 
—cAfo, where c parametrizes the decay likelihood per unit of time. This equation has the solution A^(t) = 
A/o(l — ct), which becomes negative for i > 1/c, because by then the probability for having had a 
decay exceeds unity. So what we made wrong was not to take into account that only an undecayed 
nucleus can decay, i.e. that the equation ought to have been dM /dt = —cM (t) with the solution M{t) = 
Aoexp(— cf). This is a nicely well-behaved expression, where the total probability for decays goes to 
unity only for t ^ oo. If c had not been a constant but varied in time, c = c{t), it is simple to show that 
the solution instead would have become 



Mo exp c{t') dt^ ^ ^ = -c{t)NQ exp ^- ^ c{t') dt^ . (17) 



For a shower the relevant "time" scale is something like 1/Q, by the Heisenberg uncertainty 
principle. That is, instead of evolving to later and later times we evolve to smaller and smaller Q"^. 
Thereby the DGLAP eq. (10) becomes 

^^"-^^ (- ^ C^^ ""-'^^''^ ^^') ' ^''^ 

where the exponent (or simple variants thereof) is the Sudakov factor. As for the radioactive-decay 
example above, the inclusion of a Sudakov ensures that the total probability for a parton to branch never 
exceeds unity. Then you may have sequential radioactive decay chains, and you may have sequential 
parton branchings, but that is another story. 

It is a bit deeper than that, however. Just as the standard branching expressions can be viewed 
as approximations to the complete matrix elements for real emission, the Sudakov is an approximation 
to the complete virtual corrections from loop graphs. The divergences in real and virtual emissions, so 
strange-looking in the matrix-element language, here naturally combine to provide a physical answer 
everywhere. What is not described in the shower, of course, is the non-universal finite parts of the real 
and virtual matrix elements. 

The implementation of a cascade evolution now makes sense. Starting from a simple qq system 
the q and q are individually evolved downwards from some initial Q^ax until they branch. At a branching 
the mother parton disappears and is replaced by two daughter partons, which in their turn are evolved 
downwards in and may branch. Thereby the number of partons increases, until the lower cutoff scale 
is reached. 

This does not mean that everything is uniquely specified. In particular, the choice of evolving in 
q2 _ ^2 means obvious. Any alternative variable = f{z) Q"^ would work equally well, 

since dP'^/P^ = dQ'^/Q'^. Alternative evolution variables therefore include the transverse momentum, 

z{l — z)m?, and the energy-weighted emission angle £"^6"^ m? / (z(l — z)). 
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Both these two alternative choices are favourable when the issue of coherence is introduced. Co- 
herence means that emissions do not occur independently. For instance, consider g\ g2 gs, followed 
by an emission of a gluon either from 2 or 3. When this gluon is soft it caimot resolve the individual 
colour charges of §2 and (73, but only the net charge of the two, which of course is the charge of gi. 
Thereby the multiplication of partons in a shower is reduced relative to naive expectations. As it turns 
out, evolution in p± or angle automatically includes this reduction, while one in mass does not. 

In the study of FSR, e.g. at LEP, three algorithms have been commonly used. The Herwig 
angular-ordered and Pythia mass-ordered ones are conventional parton showers as described above, 
while the Ariadne [19] p^-ordered one is based on a picture of dipole emissions. That is, instead of 
considering a — > 6 c one studies ab cde. One aspect of this is that, in addition to the branching parton, 
Ariadne also expUcitly includes a "recoil parton" needed for overall energy-momentum conservation. 
Additionally emissions off a and b are combined in a well-defined manner. 

All three approaches have advantages and disadvantages. As already mentioned, PYTHIA does 
not inherently include coherence, but has to add that approximately by brute force. Both PYTHIA and 
Herwig break Lorentz invariance slightly. The Herwig algorithm cannot cover the full phase space 
with it emissions, but has to fill in some "dead zones" using higher-order matrix elements. The ARIADNE 
dipole picture does not include g ^ qq branchings in a natural way. And so on. 

When all is said and done, it turns out that all three algorithms do a quite decent job of describing 
LEP data [20], but typically ARIADNE does best and HERWIG worst. Since Ariadne uses Pythia for 
hadronization the difference between those two is entirely due to the shower algorithms, while compar- 
isons with Herwig also are compUcated by significant differences in hadronization. 

4.3 Initial-state radiation 

The structure of initial-state radiation (ISR) is more complicated than that of FSR, since the nontrivial 
structure of the incoming hadrons enter the game. A proton is made up out of three quarks, uud, plus 
the gluons that bind them together. This picture is not static, however: gluons are continuously emitted 
and absorbed by the quarks, and each gluon may in its turn temporarily split into two gluons or into a qq 
pair. Thus a proton is teeming with activity, and much of it in a nonperturbative region where we cannot 
calculate. We are therefore forced to introduce the concept of a parton density fb{x, Q'^) as an empirical 
distribution, describing the probability to find a parton of species 6 in a hadron, with a fraction x of the 
hadron energy-momentum when the hadron is probed at a resolution scale Q^. 

While fb{x, Q'^) itself cannot be predicted, the change of with resolution scale can, once is 
large enough that perturbation theory should be apphcable: 



This is actually nothing but our familiar DGLAP equations. Before they were written in an exclusive 
manner: given a parton a, what is the probability that it will branch to 6 c during a change dQ^l Here 
the formulation is instead inclusive: given that the probability distributions fa(x, Q^) of all partons a are 
known at a scale Q^, how is the distribution of partons b changed by the set of possible branchings a b 
(+c, here implicit). The splitting kernels Pa^bci^) are the same to leading order, but differ between ISR 
and FSR in higher orders. In higher orders also the concept of fb{x, Q^) as a positive definite probability 
is lost, additional complications that we will not consider any further here. 

Even though eqs. (10) and (19) are equivalent, the physics context is different. In FSR the outgoing 
partons have been kicked to large timelike virtualities by the hard process and then cascade downwards 
towards the mass shell. In ISR we rather start out with a simple proton at early times and then allow 
more and more spacelike virtualities as we get closer to the hard interaction. Not that big fluctuations 
could not happen at early times — they do — but if they happen too early the uncertainty relation does 




(19) 
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Fig. 5: A cascade of successive branchings. The thick line represents the main chain of spacelike partons leading 
in to the hard interaction (marked by a cross). The thin lines are partons that cannot be recombined, while dashed 
lines are further fluctuations that may (if spacelike) or may not (if timelike) recombine. In this graph lines can 
represent both quarks and gluons. 

not allow them to live long enough to be of any intererest to us. That is, the higher the virtuality, the later 
the fluctuation has to occur. 

So, when the hard scattering occurs, in some sense the initial-state cascade is already there, as a 
virtual fluctuation. Had no collision occured the fluctuation would have collapsed back, but now one of 
the partons of the fluctuation is kicked out in a quite different direction and can no longer recombine with 
its sister parton from its last branching, nor with its aunt from the last-but-one branching. And so on for 
each preceding branching in the cascade that lead up to this particular parton. Post facto we therefore see 
that a chain of branchings with increasing values built up an ISR shower. Fig. 5. 

The obvious way to simulate this situation would be to pick partons in the two incoming hadrons 
from parton densities at some low scale, and then use the exclusive formulation of eq. (10) to construct 
a complete picture of partons available at higher scales, event by event. The two sets of incoming 
partons could then be weighted by the cross section for the process under study. A problem is that this 
may not be very efficient. We have to evolve for all possible fluctuations, but at best one particular parton 
will collide and most of the other fluctuations will collapse back. The cost may become prohibitive when 
the process of interest has a constrained phase space, like a light-mass Higgs which has to have the 
colliding partons matched up in a very narrow mass bin. 

There are ways to speed up this "forwards evolution" approach. However, the most common 
solution is instead to adopt a "backwards evolution" point of view [21]. Here one starts at the hard 
interaction and then tries to reconstruct what happened "before". To be more precise, the cross-section 
formula in eq. (3) already includes the summation over all possible incoming shower histories by the 
usage of Q^-dependent parton densities. Therefore what remains is to pick one exclusive shower history 
from the inclusive set that went into the Q^-evolution. To do this, recast eq. (19) as 

dn = ^A = |d(i„e^)i j:/d.^:«fl^l|ip„^,.(. = ii) . (20) 

h ^ J xfb{x,t) 2tt \ x'J 

Then we have defined a conditional probability: if parton h is present at scale Q^, what is the probability 
that it will turn out to have come from a branching a — > 6 c at some infinitesimally smaller scale? (Recall 
that the original eq. (19) was defined for increasing virtuality.) Like for FSR this expression has to be 
modified by a Sudakov factor to preserve total probability, and this factor is again the exponent of the 
real-emission expression with a negative sign, integrated over from an upper starting scale Q^ax 
down to the of the hypothetical branching. 

The approach is now clear. First a hard scattering is selected, making use of the Q^-evolved parton 
densities. Then, with the hard process as upper maximum scale, a succession of ISR branchings are 
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reconstructed at lower and lower scales, going "backwards in time" towards the early low-virtuality 
initiators of the cascades. Again some cutoff needs to be introduced when the nonperturbative regime is 
reached. 

Unfortunately the story does not end there. For FSR we discussed the need to take into account 
coherence effects and the possibility to use different variables. Such issues exist here as well, but also 
additional ones. For instance, evolution need not be strictly ordered in Q^. (Eq. (4) only gives you 
Q\ > zQl, not Ql > (5a.) Non-ordered chains in some cases can be important. Another issue is that 
there can be so many partons evolving inside a hadron that they become close-packed, which leads to 
additional recombinations. See the lectures of G. Ingelman for further details [22]. 

In summary, ISR and FSR share many aspects, but also differ. The DGLAP evolution with Su- 
dakov factors allows a simple probabilistic framework, where an initial parton undergoes successive 
branchings. For ISR the branching (usually) is in terms of higher and higher spacelike virtuaUties as the 
hard scattering approaches, while for FSR the branchings involve lower and lower timelike virtualities 
as the hard scattering recedes. In FSR both daughter partons appear on equal footing, in that both can 
be timelike and branch further. In ISR only the daughter parton on its way in to the hard scattering can 
be spacelike virtual; its sister will become part of the final state and thus has to be on mass shell, or else 
be timehke and start an FSR shower of its own. And, last but not least, the whole FSR framework is 
considerably better understood than the ISR one. 

5 Combining matrix elements and parton showers 

As we have seen, both matrix elements (ME) and parton showers (PS) have advantages and disadvan- 
tages. 

To recall, ME allow a systematic expansion in powers of Og, and thereby offer a controlled ap- 
proach towards higher precision. Calculations can be done with several (up to ~ 8) partons in the final 
state, so long as only Born-level results are asked for, and it is possible to tailor the phase-space cuts 
for these partons precisely to the experimental needs. Loop calculations are much more difficult, on the 
other hand, and the mathematically correct cancellation between real- and virtual-emission graphs in the 
soft/collinear regions is not physically sensible. Therefore ME cannot be used to explore the internal 
structure of a jet, and are difficult to match to hadronization models, which are supposed to take over in 
the very soft/collinear region. 

PS, on the other hand, clearly are approximate and do not come with a guaranteed level of precision 
for well separated jets. You caimot steer the probabihstic evolution of a shower too much, and therefore 
the efficiency for obtaining events in a specific region of phase space can be quite low. On the other 
hand, PS are universal, so for any new model you only need to provide the basic hard process and then 
PS will turn that into reasonably realistic multiparton topologies. The use of Sudakov factors ensures a 
physically sensible behaviour in the soft/collinear regions, and it is also here that the PS formaUsm is 
supposed to be as most reliable. It is therefore possible to obtain a good picture of the internal structure 
of jets, and to provide a good match to hadronization models. 

In a nutshell: ME are good for well separated jets, PS for the structure inside jets. Clearly the two 
complement each other, and a marriage is highly desirable. To do this, without doublecounting or gaps in 
the phase space coverage, is less trivial, and several alternative approaches have been developed. In the 
following we will discuss three main options: merging, vetoed parton showers and MC@NLO, roughly 
ordered in increasing complexity. Which of these to use may well depend on the task at hand. 

5.1 Merging 

The aspiration of merging is to cover the whole phase space with a smooth transition from ME to PS. 
The typical case would be a process where the lowest-order (LO) ME is known, as well as the next-to- 
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leading-order (NLO) real-emission one, say of an additional gluon. The shower should then reproduce 



W^^ = — + (21) 

a(LO) d(phasespace) 

starting from a LO topology. If the shower populates phase space according to W^^ this implies that a 
correction factor W^^^ /W^^ need to be applied. 

At first glance this does not apper to make sense: if all we do is get back W^^, then what did we 
gain? However, the trick is to recall that the PS formula comes in two parts: the real-emission answer 
and a Sudakov factor that ensures total conservation of probability. What we have called W^^ above 
should only be the real-emission part of the story. It is also this one that we know will agree with 
in the soft and collinear regions. Actually, with some moderate amount of effort it is often possible to 
ensure that W^^^/W^^ is of order unity over the whole phase space, and to adjust the showers in the 
hard region so that the ratio always is below unity, i.e. so that standard Monte Carlo rejection techniques 
can be used. What the Sudakov factor then does is introduce some ordering variable Q^, so that the 
whole phase space is covered starting from "hard" emissions and moving to "softer" ones. At the end of 
the day this leads to a distribution over phase space hke 
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W^tUQ'') = W^'^m exp ( - / ""^ W'^^iQ") dQ" . (22) 



That is, we have used the PS choice of evolution variable to provide an exponentiated version of the ME 
answer. As such it agrees with the ME answer in the hard region, where the Sudakov factor is close to 
unity, and with the PS in the soft/collinear regions, where VF^^^ « W^^. 

The method is especially convenient for resonance decays, such as e+e" '-f*/Z^ qq where 
it was first introduced [23]. In that case there is an added bonus: the full NLO answer, with virtual 
corrections included, is known to be a^^'^ = (1 -|- Os/tt). So it is trivial to use the procedure above 

and rescale everything by ( 1 + cts /tt) to obtain a complete NLO answer. Note that the difference between 
using (t(LO) or cr(NLO) in the denominator of eq. (21) only gives a difference to 0{ag), i.e. to NNLO. 

In Pythia this approach is used for essentially all resonance decays in the Standard Model and 
minimal supersymmetric extensions thereof: t bW^, — > ud, H bb, x° QQ, Q ^ Q9, 
. . . [24]. It is also used in ISR to describe qq Y/Z'^/W'^ [25], but here the NLO corrections are 
more tricky, so the cross section remains as provided by the LO number. 

Merging is also used for several processes in Herwig, such as j* /Z^ qq, t bW^ and 
qq ^* jZ^ jW^ [26]. A special problem here is that the angular-ordered algorithms, both for FSR and 
for ISR, leave some "dead zones" of hard emissions that are kinematically forbidden for the shower to 
populate. It is therefore necessary to start directly from higher-order matrix elements in these regions. A 
consistent treatment still allows a smooth joining across the boundary. 



5.2 Vetoed parton showers 

In some sense vetoed parton showers is an extension of the merging approach above. The objective is still 
to combine the real-emission behaviour of ME with the emission-ordering-variable-dependent Sudakov 
factors of PS. While the merging approach only works for combining the LO and NLO expressions, 
however, the vetoed parton showers offer a generic approach for combining several different orders. 
Therefore it is likely to be a standard tool for many studies in the future. 

To understand how the algorithm works, consider a lowest-order process such as gg — > W^. For 
each higher order one additional jet would be added to the final state, so long as only real-emission 

graphs are considered: in first order e.g. qq W^g, in second order e.g. qq W^gg, and so on. Call 
these (differential) cross sections ctq, ai, a2, It should then come as no surprise that each ai, i > 1, 
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contains soft and coUinear divergences. We therefore need to impose some set of ME phase-space cuts, 
e.g. on invariant masses of parton pairs, or on parton energies and angular separation between them. 
When these cuts are varied, so that e.g. the mass or energy thresholds are lowered towards zero, all of 
these (Tj, i> I, increase without bounds. 

However, in the ME approach without virtual corrections there is no "detailed balance", wherein 
the addition of cross section to aj+i is compensated by a depletion of ai. That is, if you have an event 
with z jets at some resolution scale, and a lowering of the minimal jet energy reveals the presence of one 
additional jet, then you should reclassify the event from being i-jet to being i + 1-jet. Add one, subtract 
one, with no net change in J2i ^i- So the trick is to use the Sudakovs of showers to ensure this detailed 
balance. Of course, in a complete description the cancellation between real and virtual corrections is not 
completely exact but leaves a finite net contribution, which is not predicted in this approach. 

A few alternative algorithms exist along these lines. All share the three first steps as follows [27]: 

1) Pick a hard process within the ME-cuts-allowed phase-space region, in proportions provided by the 
ME integrated over the respective allowed region, ctq : cti : o"2 : . . .. Use for this purpose a fix ago larger 
than the ctg values that will be used below. 

2) Reconstruct an imagined shower history that describes how the event could have evolved from the 
lowest-order process to the actual final state. That provides an ordering of emissions by whatever shower- 
evolution variable is intended. 

3) The "best-bet" choice of as scale in showers is known to be the squared transverse momentum of the 
respective branching. Therefore a factor Wq, = nbranchings('^s(^'ii)/"so)5 provides the probabiUty that 
the event should be retained. 

Now the algorithms part way. In the CKKW-L approach the subsequent steps are: 

4) Evaluate Sudakov factors for all the "propagator" lines in the shower history reconstructed in step 
2, i.e. for intermediate partons that split into further partons, and also for the evolution of the final 
partons down to the ME cuts without any further emissions. This provides an acceptance weight VFsud = 
n "propagators" Sudakov((5bcg, Qcnd) whcrc Q^^^ is the large scale where a parton is produced by a 
branching and Q^^j either is the scale at which the parton branches or the ME cuts, the case being. 

4a) In the CKKW approach [28] the Sudakovs are evaluated by analytical formulae, which is fast. 

4b) In the L approach [29] trial showers are used to evaluate Sudakovs, which is slower but allows a 

more precise modelling of kinematics and phase space than offered by the analytic expression. 

5) Now the matrix-element configuration can be evolved further, to provide additional jets below the ME 
cuts used. In order to avoid doublecounting of emissions, any branchings that might occur above the ME 
cuts must be vetoed. 

The MLM approach [27] is rather different. Here the steps instead are: 
4') Allow a complete parton shower to develop from the selected parton configuration. 
5') Cluster these partons back into a set of jets, e.g. using a cone-jet algorithm, with the same jet- 
separation criteria as used when the original parton configuration was picked. 
6') Try to match each jet to its nearest original parton. 

7') Accept the event only if the number of clustered jets agrees with the number of original partons, and 
if each original parton is sensibly matched to its jet. This would not be the case e.g. if one parton gave 
rise to two jets, or two partons to one jet, or an original b quark migrated outside of the clustered jet. 
The point of the MLM approach is that the probability of not generating any additional fatal jet activity 
during the shower evolution is provided by the Sudakovs used in the step 4' . 

5.3 MC@NLO 

MC@NLO [30] in some respects is the most ambitious approach: it aims to get not only real but also 
virtual contributions correctly included, so that cross sections are accurate to NLO, and that NLO results 
are obtained for all observables when formally expanded in powers of ccg- Thus hard emissions should 
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Fig. 6: MC@NLO applied to a production. The region between the two curves is considered as "true" Z + 1 jet 
events, with showers added to that. The rest — LO real, NLO virtual and NLO real in the shower approximation 
— are combined to cancel singularities and then showered as simple Z events. 

again be generated according to ME, while soft and coUinear ones should fall within the PS regime. 

In simplified terms, the scheme works as follows: 

1) Calculate the NLO ME corrections to an n-body process, including n + 1-body real corrections and 
n-body virtual ones. 

2) Calculate analytically how a first branching in a shower starting from a n-body topology would pop- 
ulate n + 1-body phase space, excluding the Sudakov factor. 

3) Subtract the shower expression from the n + 1 ME one to obtain the "true" n + 1 events, and consider 
the rest as belonging to the n-body class. The PS and ME expressions agree in the soft and collinear 
limits, so the singularities there cancel, leaving finite cross sections both for the n- and n + 1-body event 
classes. 

4) Now add showers to both kinds of events. 

A toy example, for the case of Z^ production, is shown in Fig. 6. Several more complicated 
processes have been considered, such as bb, it and W^W~ production. A technical problem is that, 
although ME and PS converge in the collinear region, it is not guaranteed that ME is everywhere above 
PS. This is solved by having a small fraction of events with negative weights. 

In summary, MC@NLO is superiour in that it does provide the total cross section for a process 
to NLO accuracy, and so is essential for a set of precision tests. The real-emission n + 1-body part is 
the same as used for merging, however, so for normalized event shapes the merging approach is as valid. 
(Any differences are of higher order.) Finally, if multijet topologies need to be studied, where several 
orders may contribute, vetoed showers are more appropriate. Each tool to its task. 

6 Multiple interactions 

The cross section for 2 — > 2 QCD parton processes is dominated by t-channel gluon exchange, as we 
already mentioned, and thus diverges like dp\/p\_ for —>■ 0. Introduce a lower cut p±mm and integrate 
the interaction cross section above this, properly convoluted with parton densities. At LHC energies this 
fintCpxmin) rcachcs around 100 mb for p_Lmm = 5 GeV, and 1000 mb at around 2 GeV. Since each 
interaction gives two jets to lowest order, the jet cross section is twice as big. This should be compared 
with an expected total cross section of the order of 100 mb. (QCD is a confining theory and thus provides 



a finite total cross section, unlike QED, where infinitely small scattering angles are allowed at infinitely 
large distances.) In addition, at least a third of the total cross section is related to elastic scattering 
PP PP low-mass diffractive states pp pX that could not contain jets. 

So can it really make sense that o'int(p_Lmin) > ctot? Yes, it can! The point is that each incoming 
hadron is a bunch of partons. You can have several (more or less) independent parton-parton interactions 
when these two bunches pass through each other. And the point is that an event with n interactions above 
P±min counts oncc for the total cross section but once for each interaction when the interaction rate is 
calculated. That is, 

oo oo 

0-tot = X] ^" "^^^^^ '^i"* = ^'rian , (23) 

n=0 n=0 

where (j„ is the cross section for events with n interactions. Thus dint > ctot is equivalent to (n) > 1, 
that each event on the average contains more than one interaction. Furthermore, if interactions really 
do occur independently when two hadron pass by each other, then one would expect a Poissonian dis- 
tribution, Vn = {n}"' exp(— (n))/n!, so that several interactions could occur occasionally also when 
o'int(p±min) < o"tot» e.g. for a larger p±mm cut. Energy-momentum conservation ensures that interac- 
tions never are truly independent, and also other effects enter (see below), but the Poissonian ansatz is 
still a useful starting point. 

Multiple interactions (MI) can only be half the solution, however. The divergence for p^min 
would seem to imply an infinite average number of interactions. But what one should realize is that, in 
order to calculate the do- / di matrix elements within standard perturbation theory, it has to be assumed 
that free quark and gluon states exist at negative and positive infinity. That is, the confinement of colour 
into hadrons of finite size has not been taken into account. So obviously perturbation theory has to have 
broken down by 

h 0.2GeV-fm , 

P±min ^ - ~ — ~ 0.3 GeV ~ Aqcd . (24) 

rp 0.7 im 

The nature of the breakdown is also easy to understand: a small-px gluon, to be exchanged between 
the two incoming hadrons, has a large transverse wavelength and thus almost the same phase across the 
extent of each hadron. The contributions from all the colour charges in a hadron thus add coherently, and 
that means that they add to zero since the hadron is a colour singlet. 

What is then the typical scale of such colour screening effects, i.e. at what p± has the interaction 
rate dropped to ~half of what it would have been if the quarks and gluons of a proton had all been free to 
interact fully independently? That ought to be related to the typical separation distance between a given 
colour and its opposite anticolour. When a proton contains many partons this characteristic screening 
distance can well be much smaller than the proton radius. Empirically we need to introduce a p±min 
scale of the order of 2 GeV to describe Tevatron data, i.e. of the order of 0.1 fm separation. It is not 
meaningful to take this number too seriously without a detailed model of the space-time structure of a 
hadron, however. 

The 2 GeV number is very indirect and does not really tell exactly how the dampening occurs. One 
can use a simple recipe, with a step-function cut at this scale, or a physically more reasonable dampening 
by a factor p^/ (p^Q +P±)'^, plus a corresponding shift of the as argument. 



dp]_ p\_ {pIo+pD^ 



with p±Q a dampening scale that also lands at around 2 GeV. This translates into a typical number of 2-3 
interactions per event at the Tevatron and 4-5 at LHC. For events with jets or other hard processes the 
average number is hkely to be higher. 
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6.1 Multiple-interactions models 

The first studies of complete events based on perturbatively generated MI [31] started out from a mimi- 
mally simple model: 

1) Use a sharp p^min cut as the key tunable parameter. 

2) Only address inelastic nondiffractive events, i.e. with ~ (1/2 — 2/3)crtot> so that the average 
number of interactions per such events is (n) = crint(p±min)/cnd- 

3) To first approximation this gives a Poissonian distribution in the number of interactions per event, with 
a fraction Vq = e~^"^ of purely lov/-p± interactions. 

4) The interactions are generated in an ordered sequence of decreasing p± values: p±i > p±2 > p±3 > 
This is possible with the standard Sudakov kind of trick: 

dp 1 do- 

with a starting p±o = £'cm/2. 

5) The ordering of emissions allows parton densities to be rescaled in x after each interaction, so that 
energy-momentum is not violated. Thereby the actual distribution of the number of interactions becomes 
narrower than Poissonian, since the rate of a further interaction is reduced if previous ones already took 
away energy. 

6) For technical reasons the model was simplified after the first interaction, so that there only gg or qq 
outgoing pairs were allowed, and no showers were added to these additional 2^2 interactions. 
Already this simple PYTHlA-based model is able to "explain" a large set of experimental data. 

Some additional features were then included in a more sophisticated variant. 

1) Use the smooth turnoff of eq. (25) all the way down to p± = 0. Then an event has to contain at least 
one interaction to be an event at all, i.e. p±o has to be selected sufficiently small that a-mt > <^nd- 

2) Hadrons are extended, and therefore partons are distributed in (transverse) coordinates. To allow a 
flexible parametrization and yet have an easy-to- work- with expression, a double Gaussian Pmatter(r) = 
A^i exp (— r^/r^) + N2 exp (— r^/rl) is used, where N2/N1 and r2/r\ are tunable parameters. 

3) The matter overlap during a collision, calculated by 

0{h) = I d'^dt p\X^A^,t)p^2X^ri^,t) , (27) 

directly determines the average activity in events at different impact parameter b: {n{h)) cc 0{b). That 
is, central collisions tend to have more activity, peripheral less, but of course properly normahzed so that 
the 6-integrated interaction cross section agrees with standard perturbation theory (modulo the already- 
discussed dampening at small p±). Thereby the Vn distribution becomes broader than a Poissonian. 

As before, several simplifications are necessary. 

This is the scenario that has been used in many of the experimental studies over the years. 

More recently a number of improvements have been included [32]. 

1) The introduction of junction fragmentation, wherein the confinement field between the three quarks 
in a baryon is described as a Y-shaped topology, now allows the handUng of topologies where several 
valence quarks are kicked out, thus allowing arbitrary flavours and showering in all interactions in an 
event. 

2) Parton densities are not only rescaled for energy-momentum conservation, but also to take into ac- 
count the number of remaining valence quarks, or that sea quarks have to occur in qq pairs. 

3) The introduction of -ordered showers allows the selection of new ISR branchings and new interac- 
tions to be interleaved in one common sequence of falling p± values. Thereby the competition between 
these two components, which both remove energy from the incoming beams, is modelled more reahsti- 
cally. FSR is not yet interleaved, but also does not compete for beam energy. 

This scenario is not yet as well studied experimentally. 
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The traditional Herwig soft underlying event (SUE) approach to this issue has its origin in the 
UA5 Monte Carlo. In it a number of clusters are distributed almost independently in rapidity and trans- 
verse momentum, but shifted so that energy-momentum is conserved, and the clusters then decay isotrop- 
ically. The multiplicity distribution of clusters and their y and p_\_ spectra are tuned to give the observed 
inclusive hadron spectra. No jets are produced in this approach. 

The Jimmy [33] program is an add-on to Herwig. It replaces the SUE model with a Ml-based one 
more similar to the Pythia ones above, e.g. with an impact-parameter-based picture for the multiple- 
interactions rate. Technical differences exist, e.g. JiMMY interactions are not picked to be p_L-ordered 
and thus energy-momentum issues are handled differently. 

The DPMjet/DTUjet/PhoJet family of programs [34] come from the "historical" tradition of 
soft physics, wherein multiple ~ "pomeron" exchanges fill a role somewhat similar to the hard 
MI above. Jet physics was originally not included, but later both hard and soft interactions have been 
allowed. One strong point is that this framework also allows diffractive events to be included as part of 
the same basic machinery. 

6.2 Multiple-interactions studies 

How do we know that MI exist? The key problem is that it is not possible to identify jets coming from 
2 GeV partons. Therefore we either have to use indirect signals for the presence of interactions at 
this scale or we have to content ourselves with studying the small fraction of events where two interac- 
tions occur at visibly large values. 

An example of the former is the total charged multiplicity distribution in high-energy pp/pp col- 
hsion. This distribution is very broad, and is even getting broader with increasing energy, meaured in 
terms of the width over the average, a {rich) / (rich)- By contrast, recall that for a Poissonian this quantity 
scales like 1 / ^Jn^^ and thus is getting narrower. Simple models, with at most one interaction and with a 
fragmentation framework in agreement with LEP data, cannot explain this: they are way to narrow, and 
have the wrong energy behaviour. If MI are included the additional variability in the number of inter- 
actions per event offers the missing piece [31]. The variable-impact-parameter improves the description 
further. 

Another related example is forward-backward correlations. Consider the charged multiplicity nj 
and n(, in a forward and a backward rapidity bin, each of width one unit, separated by a central rapidity 
gap of size Ay. It is not uimatural that Uf and n?, are somewhat correlated in two-jet events, and for 
small Ay one may also be sensitive to the tails of jets. But the correlation coefficient, although falling 
with Ay, still is appreciable even out to Ay = 5, and here again traditional one-interaction models come 
nowhere near. In a MI scenario each interaction provides additional particle production over a large 
rapidity range, and this additional number-of-MI variability leads to good agreement with data. 

Direct evidence comes from the study of four-jet events. These can be caused by two separate 
interactions, but also by a single one where higher orders (call it ME or PS) has allowed two additional 
branchings in a basic two-jet topology. Fortunately the kinematics should be different. Assume the four 
jets are ordered in p±, pj^i > pj_2 > P±3 > P±4- If coming from two separate interactions the jets should 
pair up into two separately balancing sets, |p_Li -I- p_L2| ~ and |p_L3 -I- p_L4| ^ 0. If an azimuthal angle 
(f is introduced between the two jet axes this also should be flat if the interactions are uncorrelated. By 
contast the higher-order graph offers no reason why the jets should occur in balanced pairs, and the if 
distribution ought to be peaked at small values, corresponding to the familiar coUinear singularity. The 
first to observe an MI signal this way was the AFS collaboration [35] at ISR (pp at 62 GeV), but with 
large uncertainties. A more convincing study was made by CDF [36], who obtained a clear signal in 
a sample with three jets plus a photon. In fact the deduced rate was almost a factor of three higher 
than naive expectations, but quite in agreement with the impact-parameter-dependent picture, wherein 
correlations of this kind are enhanced. 
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A topic that has quite extensively studied in CDF is that of the jet pedestal [37], i.e. the increased 
activity seen in events with a jet, even away from the jet itself, and away from the recoiling jet that should 
be there. Some effects come from the showering activity, i.e. the presence of additional softer jets, but 
much of it rather finds its explanation in MI, as a kind of "trigger bias" effect, as follows. (1) Central 
colhsions tend to produce many interactions, peripheral ones few. (2) If an event has n interactions there 
are n chances that one of them is hard. Combine the two and one concludes that events with hard jets are 
biased towards central colhsions and many additional interactions. The rise of the pedestal with triggger- 
jet energy saturates once aint{p±mm = P±jet) <7nd» however, because by then events are already 
maximally biased towards small impact parameter. And this is indeed what is observed in the data: a 
rapid rise of the pedestal up to p±jet ~ 10 GeV, and then a slower increase that is mainly explained by 
showering contributions. 

In more detailed studies of this kind of pedestal effects there are also some indications of a jet 
substructure in the pedestal, i.e. that indeed the pedestal is associated with the production of additional 

(soft) jet pairs. 

In spite of many qualitative successes, and even some quantitative ones, one should not be lead to 
beheve that all is understood. Possibly the most troublesome issue is how colours are hooked up between 
all the outgoing partons that come from several different interactions. A first, already difficult, question 
is how colours are correlated between all the partons that are taken out from an incoming hadron. These 
colours are then mixed up by the respective scattering, in principle (approximately) calculable. But, 
finally, all the outgoing partons will radiate further and overlap with each other on the way out, and how 
much that may mess up colours is an open question. 

A sensitive quantity is (p_L)(nch), i-e- how the average transverse momentum of charged particles 
varies as a function of their multiplicity. If interactions are uncorrelated in colour this curve tends to 
be flat: each further interaction adds about as much p± as rich- If colours somehow would rearrange 
themselves, so that the confinement colour fields would not have to run criss-cross in the event, then the 
multiplicity would not rise as fast for each further interaction, and so a positive slope would result. The 
embarrassing part is that the CDF tunes tend to come up with values that are about 90% on the way to 
being maximally rearranged [37], which is way more than one would have guessed. Obviously further 
modelling and tests are necessary here [38]. 

Another issue is whether the p±o regularization scale should be energy-dependent. In olden days 
there was no need for this, but it became necessary when HERA data showed that parton densities rise 
faster at small x values than had commonly been assumed. This means that the partons become more 
close-packed and the colour screening increases faster with increasing colhsion energy. Therefore an 
energy-dependent p±o is not unreasonable, but also cannot be predicted. Currently the default Pythia 
ansatz is p±o{Ecm) = (2.0 GcV) {Ecm/l-8 TcV)"-^^, i.e a predicted p±Q = 2.8 GeV at 14 TeV. This 
gives a minimum-bias multiplicity of about 7 per unit of rapidity in the central region, and a pedestal 
under jet events of around 30 charged particles per unit [39]. However, these numbers are model- and 
parameter-dependent. PhoJet predicts about half as big a pedestal, and typical Jimmy tunes about 50% 
more, a priori leaving a big range of uncertainty to be resolved once LHC runs begin. 

7 Hadronization 

The physics mechanisms discussed so far are mainly being played out on the partonic level, while ex- 
perimentahsts observe hadrons. In between exists the very important hadronization phase, where all 
the outgoing partons end up confined inside hadrons of a typical 1 GeV mass scale. This phase caimot 
(so far?) be described from first principles, but has to involve some modelling. The main approaches 
in use today are string fragmentation and cluster fragmentation. These are described in the lectures of 
G. Ingelman [22], so this section will be very brief, only with a few comments. 

Hadronization models start from some ideologically motivated principles, but then have to add 
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"cookbook recipes" with free parameters to arrive at a complete picture of all the nitty gritty details. 
This should come as no surprise, given that there are hundereds of known hadron species to take into ac- 
count, each with its mass, width, wavefunction, coupUngs, decay patterns and other properties that could 
influence the structure of the observable hadronic state, and with many of those properties being poorly 
or not at all known, either from theory (lattice QCD) or from experiment. In that sense, it is sometimes 
more surprising that models can work as well as they do than that they fail to describe everything. 

The simpler initial state at an e+e~ collider, such as LEP, implies that this is the logical place 
to tune the hadronization framework to data [20], and thereafter those tunes can be appUed to other 
studies. One such is the internal structure of jets in hadron colhder, where the pattern in many respects 
is surprisingly well described. 

On the other hand, at the HERA e^p collider it has been observed that the relative amount of 
strange-particle production is only 2/3 of that at LEP, and of (anti)baryons only 1/2. This has no simple 
explanation within the string fragmentation model, so it acts as a useful reminder that we still do not 
know as much as we should. Also other examples could be provided. 

8 Summary and outlook 

In these lectures we have followed the flow of generators roughly "inwards out", i.e. from short-distance 
processes to long-distance ones. At the core Ues the hard process, described by matrix elements. It 
is surrounded by initial- and final-state showers, that should be properly matched to the hard process. 
Multiple parton-parton interactions can occur, and the colour flow is tied up with the structure of beam 
remnants. At longer timescales the partons turn into hadrons, many of which are unstable and decay 
further. This basic pattern is Ukely to remain in the future, but many aspects will change. 

One such aspect, that stands a bit apart, is that of languages. The traditional event generators, Uke 
Pythia and Herwig, have been developed in Fortran — up until the end of the LEP era this was the 
main language in high-energy physics. But now the experimental community has largely switched to 
C-i-i- and decisions have been taken e.g. at CERN to discontinue Fortran altogether. The older generators 
are still being used, hidden under C-i-i- wrappers, but this can only be a temporary solution, for several 
reasons. One is that younger experimentalists often need to look into the code of generators and tailor 
some parts to specific needs of theirs, and if then the code is in an unknown language this will not work. 
Another is that theory students who apply for non-academic positions are much better off if their resumes 
say "expert on object-oriented programming" rather than "Fortran fan". 

A conversion program thus has begun on many fronts. Sherpa, as the youngest of the general- 
purpose generators, was conceived from the onset as a C-i~i- package and thus is some steps ahead of 
the other programs in this respect. Herwig-h- [40] is a complete reimplementation of the HERWIG, as 
is Pythia 8 of the current Pythia 6. Both conversions have taken longer than originally hoped, but 
progress is being made and first versions exist. ThePEG [41] is a generic toolkit for event generators, 
used by Herwig-i-i- and the upcoming new Ariadne. 

There are also other aspects where we have seen progress in recent years and can hope for more: 

- Faster, better and more user-friendly general-purpose matrix-element generators with an improved 

sampling of phase space. 

- New ready-made libraries of physics processes, in particular with full NLO corrections included. 

- More precise parton showers. 

- Better matching between matrix elements and parton showers. 

- Improved models for minimum-bias physics and underlying events. 

- Some upgrades of hadronization models and decay descriptions. 

In general one would say that generators are getting better all the time, but at the same time the experi- 
mental demands are also getting higher, so it is a tight race. However, given that typical hadronic final 
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states at LHC will contain hundreds of particles and quite complex patterns buried in that, it is difficult 
to see that there are any alternatives. 

As the same time as you need to use generators, you should remain critical and be on the lookout 
for bugs and bad modelUng. Bjorken already many years ago worried about the passive attitude many 
experimentalists have towards the output of generators; they "carry the authority of data itself. They look 
like data and feel hke data, and if one is not careful they are accepted as if they were data." [42]. Don't 
you fall into that trap! 
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